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Abstract 

We present studies for optimizing the next generation of ground-based imaging atmospheric Cherenkov telescopes 
(IACTs). Results focus on mid-sized telescopes (MSTs) for CTA, detecting very high energy gamma rays in the en¬ 
ergy range from a few hundred GeV to a few tens of Te V. We describe a novel, flexible detector Monte Carlo package, 
FAST (FAst Simulation for imaging air cherenkov Telescopes), that we use to simulate different array and telescope 
designs. The simulation is somewhat simplified to allow for efficient exploration over a large telescope design pa¬ 
rameter space. We investigate a wide range of telescope performance parameters including optical resolution, camera 
pixel size, and light collection area. In order to ensure a comparison of the arrays at their maximum sensitivity, we 
analyze the simulations with the most sensitive techniques used in the field, such as maximum likelihood template 
reconstruction and boosted decision trees for background rejection. Choosing telescope design parameters representa¬ 
tive of the proposed Davies-Cotton (DC) and Schwarzchild-Couder (SC) MST designs, we compare the performance 
of the arrays by examining the gamma-ray angular resolution and differential point-source sensitivity. We further 
investigate the array performance under a wide range of conditions, determining the impact of the number of tele¬ 
scopes, telescope separation, night sky background, and geomagnetic field. We find a 30-40% improvement in the 
gamma-ray angular resolution at all energies when comparing arrays with an equal number of SC and DC telescopes, 
significantly enhancing point-source sensitivity in the MST energy range. We attribute the increase in point-source 
sensitivity to the improved optical point-spread function and smaller pixel size of the SC telescope design. 
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1. Introduction 

The ground-based imaging atmospheric Cherenkov 
telescope (IACT) technique has led to significant 
progress in the field of very high energy (VFLE; E > 
100 GeV) gamma-ray astronomy over the last 25 years. 
To date, 145 sources have been detected at VHE with 
~60 sources discovered only in the last five year^j] 
IACTs allow us to study a wide range of scientific top¬ 
ics, many uniquely accessible by VHE astronomy. Cur¬ 
rent and future generations of IACTs aim to probe the 
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origins and acceleration processes of cosmic rays Q] J2j 
0 and explore the nature of black holes and their rel¬ 
ativistic jets. Other key objectives include the search 
for dark matter, axion-like particles 0 0, and Lorentz 
invariance violation j6|. This will require extensive ob¬ 
servations on a number of source classes such as pulsars 
and pulsar wind nebulae J7J, galactic binaries |8], super¬ 
nova remnants 0, active galactic nuclei □aim, and 
gamma-ray bursts Q20GO. The extragalactic sources 
can be used as “backlights” to study the attenuation on 
the extragalactic background light, useful for constrain¬ 
ing star formation history and other cosmological pa¬ 
rameters such as the Hubble constant m. 

VHE gamma rays entering the Earth’s atmosphere 
undergo e + e~ pair production, initiating electromag¬ 
netic cascades. The relativistic charged particles in the 
shower emit Cherenkov ultraviolet and optical radia- 
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tion, which is detectable at ground level. The major¬ 
ity of the emitted Cherenkov light is narrowly beamed 
along the trajectory of the gamma-ray primary in a cone 
with an opening angle of ~ 1.0°. Due to the beaming 
effect, the majority of the Cherenkov light falls within 
a Cherenkov light pool with a diameter of 200-300 m 
and a nearly constant light density. By imaging the 
Cherenkov light emitted by the shower particles, IACTs 
are able to reconstruct the direction and energy of the 
original gamma ray and to distinguish gamma rays from 
the much more prevalent cosmic-ray background. High 
resolution imaging of the Cherenkov shower offers sig¬ 
nificant benefits for IACTs by enabling a more accurate 
measurement of the shower axis which has an intrinsic 
transverse angular size of only a few arcminutes. How¬ 
ever the finite shower width and stochastic fluctuations 
in the shower development fundamentally limit the per¬ 
formance of IACTs. 

The designs of IACTs are governed by a few key fac¬ 
tors. At low energy, the number of Cherenkov photons 
compared to the night sky background necessitates a 
large <9(10-20 m) mirror diameter and high quantum 
efficiency camera. The camera must also be able to 
capture the signal very quickly since the duration of 
a Cherenkov pulse can be as short as a few nanosec¬ 
onds. The optical point-spread function (PSF) and cam¬ 
era pixel size should ideally be suitably smaller than 
the angular dimension of the gamma-ray shower. How¬ 
ever the high cost-per-pixel of camera designs used in 
current generation IACTs has generally dictated pixel 
sizes that are significantly larger (0.1°-0.2°) than the 
angular size of shower structure. Multiple viewing an¬ 
gles of the same shower offered by an array of tele¬ 
scopes drastically improves the reconstruction perfor¬ 
mance and background rejection. Finally, at high en¬ 
ergy, the sensitivity of IACTs is limited by gamma-ray 
signal statistics, requiring an array with a large effective 
gamma-ray collection area. 

The current generation of IACTs all have single-dish 
optical systems. These have small spherical mirror 
facets attached to either a spherical dish (i.e. Davies- 
Cotton (dc) iBirayn or a parabolic dish. The parabolic 
dish reduces the time spread of the Cherenkov signal 
but introduces a larger off-axis optical PSF. An interme¬ 
diate design with a spherical dish but a larger radius of 
curvature (intermediate-DC) can be used to achieve an 
improved time spread while maintaining off-axis perfor¬ 
mance CUED- These single-dish designs are appealing 
because they are relatively inexpensive, mirror align¬ 
ment is straightforward, and the optical PSF at large 
field angles is better than that of monolithic spherical 
or parabolic reflectors ED. 


The possibility of improving the PSF (especially 
off axis) and reducing the plate scale of IACTs has 
driven the study of Schwarzschild-Couder (SC) apla- 
natic telescopes with two aspheric mirror surface^] The 
improved PSF across the field of view (FoV) allows 
for more accurate surveying and mapping of extended 
sources. The reduced plate scale is highly compatible 
with new camera technologies such as Silicon photo¬ 
multipliers or multi-anode photomultiplier tubes. These 
technologies allow for a cost-effective, finely-pixelated 
image over a large FoV. Studies have been performed 
providing solutions for mirror surfaces optimized to cor¬ 
rect spherical and coma aberrations. These solutions 
are also isochronous, allowing for a short trigger coinci¬ 
dence window EDI . The first SC prototype is still being 
developed ED and has several challenges to overcome. 
In particular, the tolerances of the mechanical structure 
in the camera and mirror alignment system are relatively 
stringent, which translates to a higher cost. To provide 
comparisons at a fixed cost, our SC simulations use a 
smaller mirror area than that of the baseline DC design. 

The Cherenkov Telescope Array (CTA) is an exam¬ 
ple of a next-generation IACT observatory. CTA aims to 
surpass the current IACT systems such as H.E.S.S. Il22l . 
MAGIC |2l and VERITAS ED by an order of mag¬ 
nitude in sensitivity and enlarge the observable energy 
range from a few tens of GeV to beyond one hundred 
TeV H25I . To achieve this broad energy range and high 
sensitivity, CTA will incorporate telescopes of three dif¬ 
ferent sizes spread out over an area of ~3 km 2 . Tele¬ 
scopes are denoted by their mirror diameter as large- 
size telescopes (LSTs, ~24 m), medium-size telescopes 
(MSTs, ~12 m), and small-size telescopes (SSTs, ~4 m 
in the current design). The baseline designs for the LST 
and MST both feature a single reflector based on the 
DC optical design. Telescope designs based on dual¬ 
reflector SC optics are also being developed for both 
medium- and small-sized telescopes. The medium-size 
SC telescope (SCT) would fill a similar role to the MST 
and predominantly contribute to the sensitivity of CTA 
in the energy range between 100 GeV and 1 TeV. In this 
paper we explore a range of telescope models but focus 
primarly on the comparison of designs with character¬ 
istics similar to the MST and SCT. In the subsequent 
discussion we use MST to refer to all telescope designs 
with a primary mirror diameter of 9-12 m. DC-MST 
and SC-MST are used to specifically refer to telescopes 
with the imaging characteristics similar to the MST and 
SCT designs, respectively. 


2 Though segmented, the minor surfaces are often referred to as a 
singular mirror for brevity. 
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The baseline design of CTA includes -four LSTs, 
~30 MSTs, and ~50 SSTs. The sensitivity could be 
improved by a factor of 2-3 in the core energy range 
by expanding the MST array with an additional 24-36 
SCTs. With these additional telescopes, the combined 
MST and SCT array enters a new regime where the in¬ 
ternal effective area is comparable to the effective area 
of events landing outside the array. These so-called con¬ 
tained events have much improved angular and energy 
resolution as well as background rejection. Extensive 
work is underway to optimize the design of CTA for the 
wide range of science goals csi- The scope of previous 
studies has been primarily on a straightforward expan¬ 
sion of existing telescope designs to larger arrays. 

In this paper, we describe a novel, flexible Monte 
Carlo simulation and analysis chain. We use them to 
evaluate the performance of CTA-like arrays over a 
large range of telescope configurations and design pa¬ 
rameters. Section 2 describes this simulation and the 
simplified detector model. In Section 3, we explain the 
analysis chain, including a maximum likelihood shower 
reconstruction using simulated templates. This recon¬ 
struction was used for comparisons between the maxi¬ 
mum sensitivity for each array configuration. In Section 
4, we show comparisons between possible CTA designs, 
focusing primarily on the number of telescopes and the 
DC versus SC designs. We conclude in Section 5. 


2. Simulation 

We have studied the performance of a variety of array 
geometries and telescope configurations for a hypothet¬ 
ical CTA site at an altitude of 2000 m. Details of the 
site model and array geometry are described in Sections 


performed using a simplified detector model described 
in Section [231 

2.1. Air-Shower Simulations 

Simulations of the gamma-ray and cosmic-ray air 
shower cascades were performed with the CORSIKA 
v6.99 Monte Carlo (MC) package |[26l and the QGSJet 
11-03 hadronic interaction model ea. We used a site 
model with an elevation of 2000 m, a tropical atmo¬ 
spheric profile, and an equatorial geomagnetic field con¬ 
figuration with (B x , B z ) = (27.5 //f, - 15.0 pT). This 
site model is identical to the one used in m and has 
similar characteristics to the southern hemisphere sites 
proposed for CTA. 

Gamma-ray showers were simulated as coming from 
a point on the sky at 20° zenith angle and 0° azimuth 


2.1 and 2.2 Simulations of the telescope response were 


angle, as measured from the local magnetic north over 
the energy range from 10 GeV to 30 TeV. Protons and 
electrons were simulated with an isotropic distribution 
that extends to 8° and 5° respectively from the direction 
of the gamma-ray primary. We use the spectral parame- 
terizations for proton and electron fluxes from [[281 . To 
account for the contribution of heavier cosmic-ray nu¬ 
clei we increase the proton flux by a factor 1.2. 

2.2. Array Geometry 

Proposed designs for CTA employ three telescope 
types (SST, MST, and LST) with variable inter¬ 
telescope spacing from 120 m to more than 200 m ns. 
The number of telescopes of each type and their separa¬ 
tions are chosen to optimize the differential sensitivity 
over the full energy range of CTA. IIT8I found that two 
balanced arrays (arrays E and I) that have 3—4 LSTs, 
18-23 MSTs, and 30-50 SSTs of 7 m diameter provide 
the best compromise in performance over the full en¬ 
ergy range of CTA while keeping the total cost of the 
array within the projected CTA budget. 

For this study we simulated an array geometry which 
is similar to the one used for MSTs and LSTs in ar¬ 
rays E and I. The array is composed of 61 telescopes 
arranged on a grid with constant inter-telescope spac¬ 
ing of 120 m (see Figure |T]>. A telescope spacing of 
about 120 m is well motivated by the characteristic size 
of the Cherenkov light pool for gamma-ray air show¬ 
ers and guarantees that multiple telescopes will sam¬ 
ple the shower within the shower light pool. Subsets 
of telescopes from the baseline array were used to con¬ 
struct arrays with a reduced number of telescopes by 
removing successive rings of telescopes along the array 
perimeter. These reduced arrays have between 5 and 
41 telescopes and encompass arrays that are similar in 
telescope number to both current IACT arrays (N te i = 
5) and the array designs currently considered for CTA 
(N te i = 25—41). We also examined the performance 
of arrays with smaller and larger inter-telescope sepa¬ 
rations (60 m-200 m) by rescaling the inter-telescope 
separation of our baseline array. 

All simulations were performed with homogeneous 
arrays composed of a single telescope type. We primar¬ 
ily consider telescope models with mirror areas between 
the current MST and LST designs. Because our study is 
focused on the performance of arrays in the core CTA 
energy range (100 GeV - 10 TeV) we did not consider 
SSTs. 

2.3. Detector Model 

Simulations of IACT arrays have traditionally been 
performed with highly detailed detector models that use 
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Figure 1: Physical telescope positions for the five array geometries 
used for this study. All geometries are composed of telescopes ar¬ 
ranged on a uniform grid with 120 m spacing. The smallest array 
is composed of five telescopes (black circles). The larger arrays are 
constructed by the addition of successive rings of telescopes around 
the array boundary up to a maximum of 61 telescopes in the baseline 
array geometry. 


optical ray-tracing to track the trajectory and time of ar¬ 
rival of individual Cherenkov photons. Because these 
models have a very large number of parameters, a brute 
force optimization of the telescope design presents a 
significant computational challenge. In order to effi¬ 
ciently study the telescope design parameter space, we 
have developed a simplified telescope simulation tool, 
FAST (FAst Simulation for imaging air cherenkov Tele¬ 
scopes), that is not tied to any particular mirror config¬ 
uration or camera technology. In the FAST model, the 
telescope characteristics are fully described by the fol¬ 
lowing parameters: 

• Effective light collection area: A opt 

• 68% containment radius of the optical PSF: R ps f 

• Camera pixel size: D p i x 

• Effective camera trigger threshold: 7\ h 

• Single photo-electron (PE) charge resolution: <x spe 

• Pixel read-noise: crb 

• Effective integration window: AT 


While this simplified model lacks the level of detail pro¬ 
vided by other simulation tools, the performance of a 
realistic telescope design can be approximated by an 
appropriate choice of these model parameters. In this 
section we describe in detail the implementation of our 
model and how each of these parameters influence the 
telescope response. 

The geometrical model of the telescope consists of 
a primary mirror of diameter D with physical mirror 
area = ir(D/2) 2 . All Cherenkov photons that in¬ 
tersect with the primary mirror surface are propagated 
through the telescope simulation. The photons collected 
by the primary mirror are folded with a wavelength de¬ 
pendent photon detection efficiency, e(A), that models 
losses from all elements in the optical system and cam¬ 
era (mirrors, lightguides, and photosensors). Applying a 
detection probability to each collected photon, we con¬ 
struct a list of detected photoelectrons (PEs) which are 
used as input to the simulation of the telescope trigger, 
camera, and optics. 

We quantify the total light-collecting power of a tele¬ 
scope by its effective light collection area, A opt (/l) = 
Av/f(/l), the product of the physical mirror area with the 
total photon detection efficiency at wavelength A. We 
compute a wavelength-averaged effective area by fold¬ 
ing A opt (/l) with a model for the wavelength distribution 
of Cherenkov light, 

A op t = [ ' P(A,z)A opt (A)dA, (1) 

J Aq 

where 

P(A,z) °c e~ T<X ' z) A~ 2 (2) 

is the normalized wavelength distribution of Cherenkov 
light at the ground for an emission altitude z and an op¬ 
tical depth for atmospheric extinction t(A,z). We use 
an atmospheric extinction model generated with MOD- 
TRAN |29l for the tropical atmosphere and an aerosol 
layer with a visibility of 50 km. For all further evalua¬ 
tions of A opt we use z =10 km and an integration over 
wavelength from 250 nm to 700 nm. 

We define a benchmark telescope with a D - 12 m 
primary diameter and a total photon detection efficiency 
that includes losses from mirror reflections and photo¬ 
sensor efficiency. We use a photosensor model with a 
spectral response that is characteristic of photomulti¬ 
plier tubes and has a peak efficiency of 24% at 350 nm. 
Losses from mirror reflections are evaluated for a single 
optical surface using a wavelength-dependent reflectiv¬ 
ity with a peak efficiency of 89% at 320 nm. This reflec¬ 
tivity is similar to that of the aluminum and aluminized 
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Figure 2: Effective light collection area versus wavelength for the 
benchmark telescope model with A opt = 11.18m 2 . The dashed black 
line shows the spectral shape for Cherenkov light emitted at an eleva¬ 
tion of 10 km after absorption by the atmosphere. 

glass mirrors used in current generation IACTs. Figure 
[2]shows the optical effective area of the telescope model 
as a function of wavelength. The effective light collec¬ 
tion area of our benchmark telescope is 11.18m 2 which 
is representative of medium-sized IACTs with ~10 m 
aperture and 50-100 m 2 mirror area. The response of 
telescopes with larger or smaller light collection areas 
is modeled using the same spectral response and mirror 
area as the benchmark telescope model but scaling the 
photon detection efficiency by the ratio A opt /l 1.18 m 2 . 

The imaging response of the telescope optical system 
is simulated by applying a model for the optical point- 
spread-function (PSF) to the distribution of true pho¬ 
ton arrival directions in the camera image plane. Af¬ 
ter applying a survival probability for detection, each 
Cherenkov photon is assigned a random offset drawn 
from the optical PSF. We parameterize the optical PSF 
as a 2D Gaussian with a 68% containment radius, R ps f, 
that is constant across the FoV. We consider values of 
R ps f between 0.02° and 0.08° which is comparable to 
the range of PSF spot sizes for the CTA telescope de¬ 
signs at both small and large field angles. All telescopes 
are simulated with an 8 deg FoV with a light collection 
area that is constant with field angle. 

Telescopes are simulated with a camera geometry 
composed of square pixels of angular width D pix that 
uniformly tile the camera FoV. Each pixel is assigned 
a time integrated signal that is the sum of the detected 
Cherenkov photons, night-sky background (NSB) pho¬ 
tons, and detector noise. The number of NSB photons 
is drawn from a Poisson distribution where the average 


(jib) is computed using an implicit time integration win¬ 
dow (AT) of 16 ns. The mean number of NSB photons 
per pixel for a telescope with effective light collection 
area A opt and pixel solid angle AG is 

Ri b = ATAGe J F nsb (A)A opt (A)dA, (3) 

where F nsb (A) is the differential NSB flux versus wave¬ 
length. We use the NSB spectral model from If30l which 
is representative of the sky brightness of an extragalac- 
tic observation field. When folded with the optical effi¬ 
ciency of our benchmark telescope model, the integral 
flux of detected NSB photons is 365 MHz deg~ 2 m 2 . 
Our benchmark telescope model has an NSB surface 
density in the image plane (Z nsh ) of 65.4 deg~ 2 for an 
integration window of 16 ns. We model the photo¬ 
sensor single photoelectron response with a Gaussian 
with cr spe = 0.4 PE. Each channel is simulated with a 
Gaussian readout noise (<x fc j of 0.1 PE. For the range 
of pixel sizes and optical throughputs considered in this 
study, the readout noise is a subdominant component of 
the pixel noise relative to NSB and is therefore not ex¬ 
pected to have a significant impact on the telescope per¬ 
formance. Fig. [3] shows simulated camera images for 
telescope models with two different pixel sizes observ¬ 
ing the same 1 TeV gamma-ray shower. 

The trigger system of an IACT array rejects noise- 
induced events while maintaining high efficiency for 
cosmic-ray signals. We simulate a two-stage trigger sys¬ 
tem composed of a camera-level trigger for each tele¬ 
scope and an array-level trigger that combines the cam¬ 
era triggers of multiple telescopes to form the final trig¬ 
ger decision. Camera trigger designs used by current 
generation IACTs and envisioned for CTA are gener¬ 
ally based on a multi-level hierarchy whereby trigger 
information from individual pixels or camera subfields 
is combined to form the camera-level trigger decision 
OTl [32 , 25]. The rate of accidental triggers is sup¬ 
pressed by requiring a time coincidence of triggers from 
neighboring pixels or camera regions. 

A useful quantity for characterizing the performance 
of different camera trigger designs is the effective cam¬ 
era threshold, the true gamma-ray image amplitude in 
PEs at which the camera trigger is 50% efficient. To 
first order the efficiency of the camera trigger depends 
only on the surface brightness of the Cherenkov shower 
image. Because the angular size of the shower is only 
a weak function of distance and energy, we can approx¬ 
imate the response of a camera trigger by applying a 
fixed threshold on the true number of Cherenkov PEs in 
the camera FoV. 
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Figure 3: Camera images of the same 1 TeV gamma-ray shower with an impact distance of 120 m simulated with two different telescope pixel 
sizes: Z) p j x = 0.16° (left) and D p j x = 0.06° (right). Both telescope models have A opt = 11.18 nr and R ps f = 0.02°. The color scale denotes the 
measured signal amplitude in PEs for each pixel. The white cross and solid line show the direction of the gamma-ray primary and the projection of 
its trajectory to the telescope image plane, respectively. 


We simulate the camera trigger by applying a thresh¬ 
old 7* on the true number of Cherenkov PEs detected 
in the camera FoV. For showers that trigger one or more 
telescopes, the array-level trigger is simulated requiring 
a multiplicity of at least two triggered telescopes. The 
camera threshold provides a single parameter model that 
we use to explore influence of the trigger threshold on 
the array-level performance. By calibrating T t h to the 
effective camera threshold of a given trigger design, we 
can also approximate the trigger response that would be 
obtained with a more detailed trigger simulation imple¬ 
mentation. 

Studies performed with the sim_telarray detector 
simulation package j28l have shown that camera trigger 
designs currently considered for the MSTs can achieve 
effective trigger thresholds of 60-80 PE for a single tele¬ 
scope accidental trigger rate of 1-10 kHz. We adopted 
a trigger threshold of 60 PE for our baseline telescope 
model with Aopt = 11.18m 2 which is comparable to the 
effective threshold of the prod-2 MST model lf33l . To 
model the effective trigger threshold for telescopes with 
different light collection areas, we used a simple scal¬ 
ing formula that approximates the threshold needed to 
maintain a constant rate of accidental triggers. If the to¬ 
tal pixel noise is dominated by NSB photons, the rate of 
accidental triggers should be proportional to the RMS 


fluctuations in the number of NSB photons collected in 
a trigger pixel which scales as A op “ if the angular pixel 
size is held fixed. Telescopes with larger effective light 
collection area achieve a lower trigger threshold through 
the suppression of these NSB fluctuations relative to the 
signal amplitude which increases linearly with A opt . We 
assign the effective trigger threshold for a telescope with 
light collection area A opt as. 


T th = 60 PE 


A opt 

11.18 m 2 


1/2 


(4) 


For the studies presented in Section [TJ we consider 
a benchmark array (M61) with 61 identical telescopes 
with A opt = 11.18 m 2 , R ps f = 0.02°, D v ; x = 0.06°, and 
7* =60 PE. Our baseline telescope model is repre¬ 
sentative of a generic medium-sized telescope design 
with SC-like imaging characteristics. In Section |4~2| we 
additionally consider other telescope models that were 
specifically chosen to match the characteristics of the 
proposed CTA telescope designs. 


2.4. Simplifications of the Detector Simulation 

The FAST package uses a highly simplified model 
of the telescope optics and camera. This allows us to 
perform a more general exploration of the IACT de¬ 
sign parameter space without focusing on the details of 
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any specific optical design or camera technology. These 
simplifications also make the FAST simulation much 
less computationally intensive than traditional simula¬ 
tion tools such as sim_telarray which further facili¬ 
tates the exploration of a large phase space of telescope 
and array design concepts. When simulating compa¬ 
rable telescope designs with FAST and sim_telarray 
we observe an order-of-magnitude reduction in compu¬ 
tation time. Here, we discuss the limitations of the ap¬ 
proach taken in the FAST simulation package and de¬ 
scribe the areas in which a more detailed simulation of 
the telescope camera and optics could potentially affect 
our results. 

FAST does not use raytracing to account for shadow¬ 
ing of the camera by telescope structure and assumes a 
simple 2D Gaussian PSF that is constant over the FoV. 
Full raytracing simulations can be used to model effects 
such as shadowing by the telescope structure and light 
losses from gaps in the mirror surfaces. Raytracing also 
allows a more realistic modeling of the telescope PSF. 
Optical aberrations intrinsic to the design of IACTs in¬ 
troduce a strong field-angle dependence to both the size 
and shape of the PSF. In the FAST telescope model, all 
effects that influence the optical performance of the tele¬ 
scope are folded into the effective optical area (A opt ) and 
the 68% containment radius of the PSF (R ps f). Shadow¬ 
ing and other light losses in the optical system are thus 
taken into account by a reduction in the effective optical 
area. The field-angle dependence of the PSF is stud¬ 
ied by comparing the performance of telescope designs 
with the best and worst PSF at any field angle. The per¬ 
formance of a real telescope should always fall between 
that of telescopes with smallest and largest PSF in the 
FoV. We did not study the effect of an asymmetric PSF 
but it is plausible to assume that the array performance 
with telescopes that have an asymmetric PSF can be es¬ 
timated by enlarging the symmetric PSF in FAST. 

FAST does not simulate the timing of signals and as¬ 
sumes an ideal data acquisition system that is only lim¬ 
ited by the irreducible noise from NSB photons. Effects 
that distort the measurement of the pixel charge such 
as cross-talk, after-pulsing, timing jitter, non-linearity, 
and saturation might lead to a reduction in performance. 
However these effects should affect all telescope types 
in a similar way in the sense that photon charges might 
be only partially reconstructed. Per the telescope perfor¬ 
mance requirements set forward by CTA, we expect the 
influence of these effects to be sub-dominant to the irre¬ 
ducible limitations on the IACT technique set by shower 
physics (NSB and shower fluctuations). Many of these 
effects can be partially mitigated with pixel-level cali¬ 
bration or pre-processing analysis procedures. For in¬ 


stance, suppression algorithms such as clipping of pixel 
signal amplitudes or the removal of isolated high PE 
pixels have been demonstrated as efficient techniques 
to reduce the impact of after-pulsing. In principle these 
effects could also be approximated in the FAST detector 
model by increasing the pixel-level noise or worsening 
the charge resolution. We note that our analysis does 
not use shower timing parameters or shower develop¬ 
ment timing in any explicit way and thus is robust with 
respect to changes in timing requirements. 

Another simplification in FAST is the trigger thresh¬ 
old decision logic that assumes any shower image above 
a certain threshold will trigger the telescope. Realis¬ 
tic trigger electronics might have several characteristics 
that have to be simulated in detail but in the early plan¬ 
ning state of an array our method is very valuable to give 
a solid estimate of the array performance as long as the 
trigger threshold is not chosen aggressively and the re¬ 
quirements of the telescope array are met. Indeed we 
show in Section 4.3 that our telescope threshold might 
be too conservative given that more realistic simulations 
find enhanced performance near the energy threshold of 
simulated arrays. For relative comparisons of arrays our 
simplifications are unimportant. To estimate the influ¬ 
ence on the absolute performance impacts of our sim¬ 
plified simulations we compare our results to a much 


more sophisticated detector simulation in Section 4.3 


3. Analysis 

The analysis of the telescope image data is performed 
using well established techniques for the analysis of 
IACT data. The analysis is performed in three stages: 
preparation of the telescope images, reconstruction of 
the event properties, and training and optimization of 
cuts. 

3.1. Image Cleaning and Parameterization 

The image analysis is applied to the telescope pixel 
amplitudes to derive a set of telescope-level parameters 
which characterize the distribution of light in each tele¬ 
scope. Analysis of the telescope image data begins with 
the application of an image cleaning analysis that se¬ 
lects pixels that have a signal amplitude that is larger 
than noise. Traditionally image cleaning has been per¬ 
formed using variations of a nearest-neighbor algorithm 
m. a search is performed for groups of neighboring 
pixels which exceed a threshold defined in terms of the 
absolute amplitude or the amplitude relative to the RMS 
noise in the pixel. These algorithms work well as long 
as the dimension of the pixel is of the same order as the 
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Cherenkov image size. However in the limit of small 
pixel sizes these algorithms will lose efficiency for low 
energy showers where the signal is spread out over too 
many pixels to be discernible above noise when only 
considering nearest neighbors. 

In order to circumvent the limitations of the nearest- 
neighbor pixel algorithms, we use an Aperture cleaning 
algorithm that performs a smoothing over the camera 
with an angular scale (R = 0.12°) that is of the same 
order as the width of a gamma-ray induced Cherenkov 
shower (0.1-0.2°). 

In order to detect efficiently images that lie on pixel 
boundaries we divide each pixel into N x N subpixels 
where N = [7J p j x /0.06 0 ]. We compute the image inten¬ 
sity in the neighborhood of subpixel i as 


s(R) = ^ SjWij(R), (5) 

i 

where w,j(R) is the fraction of the solid angle of pixel 
j contained within the circular aperture of radius R cen¬ 
tered on subpixel i (see Figure [4|. The pixel image 
threshold is defined relative to the expected noise within 
the pixel aperture 


cr(R) = 


+ Rb)wij(R) 


1/2 


( 6 ) 


For the present analysis we adopt an image threshold 
of s/cr = 7 which for our baseline array corresponds 
to an image intensity of 319 PE deg~ 2 and an inte¬ 
grated charge of 14.4 PE within the cleaning aperture. 
Any pixel for which one or more subpixels exceeds the 
cleaning threshold is flagged as an image pixel. The 
simulations do not include photodetector after-pulsing, 
which can cause noise isolated in single pixels. These 
may need to be suppressed if the aperture cleaning 
method is applied in other scenarios. Telescope images 
are discarded at this point if fewer than three image pix¬ 
els are present. 

The image cleaning is only used by the geometric re¬ 
construction, itself a seed for the likelihood reconstruc¬ 
tion. As such, a relatively low threshold was chosen to 
maximize the reconstruction efficiency for low-energy 
events. 

Following the image cleaning analysis, an image 
analysis is applied to the amplitudes of image pixels (sj) 
to calculate a set of image parameters that characterize 
the light distribution in the focal plane. The image pa¬ 
rameters include the total image size, S , the image cen¬ 
troid, the second central moments along the major and 
minor axes of the image denoted as length / and width 
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Figure 4: Illustration of the aperture cleaning algorithm on small cam¬ 
era subsections with R = 0.12°. In the DC-like case (left), pixels are 
subdivided since they are large compared to the aperture. Each sub¬ 
pixel is used as the center of an aperture for image intensity calcula¬ 
tion. This calculation is based on the number of PEs and the fraction 
of the pixel area within the aperture, normalized to the area of the 
aperture. For the SC-like case (right), smaller pixels do not require 
subdivision. 


w, and the orientation of the major axis in the image 
plane. 

3.2. Shower Reconstruction 

The shower reconstruction determines a trajectory 
and energy for each event by fitting a shower model to 
the telescope image data. The shower model parame¬ 
ters ( 0 ) are the primary energy (£), the primary direc¬ 
tion (e), the primary impact position (R), and the atmo¬ 
spheric column depth of the first interaction point (A). 
In an array of IACTs, each telescope views the shower 
from a different perspective and provides an indepen¬ 
dent constraint on the shower parameters. By using 
image data from multiple telescopes, one can perform 
a stereoscopic reconstruction of the shower trajectory. 
For the analysis algorithms presented in this section, we 
assume on-axis observations of a gamma-ray source in 
parallel pointing mode whereby the optical axes of the 
telescopes in the array are aligned with the shower di¬ 
rection. However the procedures described here can be 
also applied to the case of non-aligned telescope point¬ 
ing. 

In presenting the implementation of the shower re¬ 
construction algorithms, we use a global coordinate sys¬ 
tem defined with the x-axis parallel to the direction 
of magnetic north and the z-axis perpendicular to the 
Earth’s surface. The positions of the array telescopes 
are denoted by r,. For the array layouts considered for 
this study, the telescopes are arranged in a regular grid 
in the x-y plane with all telescopes located at the same 
height above sea level (z = 2000 m). Shower recon¬ 
struction is performed in a shower coordinate system 
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with the z-axis aligned with the shower trajectory and 
defined by the basis vectors: 


z' = e 


/v/ 

y 


(e-(z-e)z) „ 

, - xz 

\/l - (z • e) 2 


x = i' x y 


(7) 


We use r' and R' to represent the projections of the tele¬ 
scope positions and the shower impact position to the 
x' -y' plane. An illustration of the geometry of a shower 
is shown in Figure [5] The shower impact vector. 


Pi = R-r ; -(z'-(R-r,■))!', (8) 

describes the location of the shower impact position rel¬ 
ative to telescope i in the x' - y' plane. The shower 
impact distance (p ,■ = |p ; [) is the distance of closest ap¬ 
proach between the shower and the telescope. 

The geomagnetic field (GF) alters the development 
of the gamma-ray shower by deflecting the charged 
particles in the electromagnetic cascade. The GF de¬ 
flects particles in a plane perpendicular to their trajec¬ 
tories with a strength proportional to the perpendicular 
component of the GF vector. For the shower particles 
that predominantly contribute to the emitted Cherenkov 
light, the perpendicular component is comparable to the 
GF vector component perpendicular to the shower di¬ 
rection (B ± = B - (B • z')z'). Deflection of the shower 
particles by the GF breaks the azimuthal symmetry of 
the shower causing an elongation in the shower particle 
distribution in the plane orthogonal to B ± . 

Due to the asymmetry in the shower development in¬ 
duced by the GF, the Cherenkov light distribution ob¬ 
served by a telescope depends on both the distance to 
the shower impact position (p) and the orientation of 
the shower impact vector relative to the GF. We param¬ 
eterize the shower orientation with respect to telescope 
i by the shower position angle (</>;) defined by 


cos (pj — x! ' 7 ~y • (9) 

\Pi\ 

Telescopes with a shower position angle of 0° and 90° 
view the shower in the planes parallel and perpendicular 
to its elongated axis respectively (see Figure 5). 

The shower reconstruction is performed in two con¬ 
secutive stages. A geometric reconstruction algorithm 
is first used to obtain a robust estimation of the shower 
parameters. In this stage the shower energy and inter¬ 
action depth are initially assigned using look-up tables. 
In the second stage the shower parameters derived from 



Figure 5: Illustration of the geometry of a gamma-ray shower as 
shown in the shower coordinate system. The gamma-ray trajectory 
is defined by its impact position R' in the x' - y' plane and arrival 
direction e. The GF induces an elongation in the shower in the plane 
orthogonal to Bj_ (indicated by the grey shaded square). The shower 
impact vector, p,, describes the position of the shower impact position 
relative to the telescope at r, (closed blue circle). The shower position 
angle, (pi, is defined by the angle between the shower impact vector 
and the x'-axis. 


the geometric reconstruction are refined using a likeli¬ 
hood reconstruction algorithm that performs a joint fit 
to the image intensity in all telescopes. 

3.2.1. Geometric Reconstruction 

The geometric reconstruction algorithm is a 3-D 
stereoscopic reconstruction technique based on the tra¬ 
ditional Hillas image parameterization of the shower 
images (35). The emitted Cherenkov light from a 
gamma-ray shower produces an approximately ellipti¬ 
cal distribution in the telescope focal plane with the ma¬ 
jor axis of the ellipse aligned with the shower trajec¬ 
tory. The projected shower trajectory as observed by a 
telescope with impact vector p, can be described by the 
equation 




Pi + e f 

I Pi + ef| ’ 


( 10 ) 


where t is the physical distance along the shower axis 
from its intersection with the shower plane. Each tele¬ 
scope that observes the shower constrains the trajectory 
to lie in the plane formed by the vectors z' and p,. When 
multiple telescope images are present, the intersection 
point of the projected shower axes provides a unique 
solution for both the shower direction (e) and its impact 
position in the shower plane (R')- 

The solution for the shower trajectory that best 
matches the projected shower axes of all telescopes is 
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found by minimizing a pair of ^-likc parameters that 
independently optimize the shower direction and core 
position. In the case of the shower direction we solve 
for the vector e that minimizes 

X 2 e (e) = J]K( s j,Wi, l i) A e /e) 2 , (11) 

i 

where A e , (e) is the distance of closest approach between 
the major axis of the image ellipse and the shower direc¬ 
tion projected to the image plane of telescope i, and k is 
a weighting function that controls the contribution of 
each telescope to the total sum. Images that are brighter 
and more elongated provide a better constraint on the 
shower trajectory, and therefore we use as our weight¬ 
ing function the product of the image size with square 
of the image ellipse eccentricity, 

ij - w? 

K(S l ,w l ,l i ) = S l ^- r ^. (12) 

I 

The shower core position is reconstructed by mini¬ 
mizing 


**(R) = 2 K(S <- w - 1 ‘)Arj(R) 2 , (13) 

i 

where 

Ajw(R) = |p,(R) - (p,(R) ■ e p ,j e p ,,| (14) 

is the distance of closest of approach between the im¬ 
age axis of telescope i projected to the shower plane 
(e p j) and the core location. After reconstruction of the 
shower trajectory, the shower energy is reconstructed 
using look-up tables for the shower energy as a func¬ 
tion of the image size and impact distance from the tele¬ 
scope. The shower energy estimate is calculated from 
a weighted average of telescope energy estimates given 

by 


E = 


V 1 

^ a- E (Si,Pi) 

i > 


sr' E(S j,pi ) 
V (T e(S i,PiY 


(15) 


where E(S,,pi) and cr £ (5,,p,) are functions for the ex¬ 
pectation value and standard deviation of the shower en¬ 
ergy derived from simulations. 


3.2.2. Likelihood Reconstruction 

The likelihood reconstruction performs a global fit to 
the telescope image data using a model for the expected 
pixel amplitude p{0) as a function of the shower param¬ 
eters 0. Pixel expectation values are evaluated from an 
image template model, /(e; p, 6), a probability distribu¬ 
tion function for the image intensity in the direction e 


as measured by a telescope that observes a shower with 
parameters 0 and impact vector p. More details on the 
generation of the image intensity model are presented 
in Section 3.2.3 The agreement between the telescope 
image model and the data is evaluated by means of an 
array likelihood function. Shower parameters are deter¬ 
mined by a maximization of an array likelihood func¬ 
tion. Maximization of the array likelihood as a function 
of shower fit parameters is performed using a numeri¬ 
cal non-linear optimization technique. In order to en¬ 
sure stable fit convergence, the shower parameters are 
initially seeded with a set of values derived by the geo¬ 
metric reconstruction ( 0 geo ). 

We use a formulation of the array likelihood function 
which is similar to the one presented in l36l . The ar¬ 
ray likelihood is computed from a pixel-by-pixel com¬ 
parison between the observed and predicted image in¬ 
tensities. The likelihood provides a statistical model 
for the measured pixel signal (,v) as a function of in¬ 
put models for signal and background. The measured 
pixel signal is modeled as the sum of three components: 
Cherenkov signal photons, NSB photons, and Gaussian 
noise arising from detector fluctuations. The pixel like¬ 
lihood function is 


Lpix(s\p(Q'), Pb, CT b , CT S p e ) — ^ ^ 


(p + Pbfe-^^ 


g(s, n). 


(16) 

where p is the model amplitude, p b is the NSB ampli¬ 
tude, cr b is the standard deviation of the detector noise, 
<x spe is the width of the single PE response function, and 


g(s,n) = 


cr 2 b + nv%e) 


exp 


(S - nf 


2 (' 


erf + ncr: 


spe \ 


(17) 

The model amplitude for pixel j in telescope i is cal¬ 
culated by an integration of the image template model 
over the pixel. 


Pi0)= f He;p„ff)dQ, (18) 

Jn,j 

where Q (/ is the 2-D angular integration region. 

The array likelihood is calculated from the product of 
the pixel likelihoods in all telescopes, 


L(s\p(0), p b , o- b ,cr y ) = Y\L pix (s i j\p ij (0),p h ,o- b ,o- y ), 
i*j 

( 19 ) 
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where s lt] and p,j are the signal and model amplitude of 
pixel i in telescope j. The set of pixels included in the 
computation of Equation [T9] can encompass the entire 
camera. Unlike for the geometric reconstruction tech¬ 
niques, each pixel is weighted by its expected contribu¬ 
tion to the total image intensity. Therefore the inclusion 
of pixels on the shower periphery does not significantly 
improve or degrade the reconstruction performance. Al¬ 
though the array likelihood can be calculated using all 
pixels in the camera, using a smaller number of pix¬ 
els significantly reduces the computation time needed 
for the shower likelihood optimization. In order to se¬ 
lect pixels that will provide a useful constraint on the 
shower parameters, we choose a set of pixels *P in each 
telescope that satisfies the relation 

2^/0geo) > /£/f,(0 g c°), (20) 

j^P j 

where p s is the expected image intensity in pixel j and / 
is the fraction of the total image intensity. We build the 
set *P by adding pixels in order of their expected inten¬ 
sity until the total amplitude fraction exceeds /. Having 
found that the reconstruction performance is relatively 
insensitive to / for values > 0.75, we use / = 0.75. 

An underlying assumption of the likelihood formula¬ 
tion presented here is that the shower can be treated as a 
continuous distribution of particles. However, because 
the electromagnetic cascade is a stochastic process, non- 
statistical deviations from the image model are expected 
due to fluctuations in the shower development. These 
deviations become especially important at low energies 
where the total number of shower particles is small and 
the influence of the GF becomes large. These shower 
fluctuations will tend to worsen the performance of the 
method relative to what would be expected in the case 
of purely statistical fluctuations. 

3.2.3. Image Templates 

The image model, /(e; p, 0), is the probability distri¬ 
bution function for the measured telescope image inten¬ 
sity in photoelectrons (PEs) versus direction, e. The 
model is parameterized as a function of the shower 
properties (energy and first interaction depth) and the 
impact position of the shower relative to the telescope. 
The model is generated by averaging the intensity of 
a large sample of simulated showers generated at a se¬ 
quence of fixed offsets, energies and interaction depths. 
The image templates for this study were generated 
with the CORSIKA shower simulation package and the 
detector simulation described in Section 12.31 While 
the image templates used for this study are MC-based 
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Figure 6: Image intensity templates for three different three gamma- 
ray energies (100 GeV, 1 TeV, and 10 TeV) generated for a telescope 
model with R ps f = 0.02° and A upl = 11.18m 2 . The images show 
the expectation for the measured intensity of Cherenkov light as a 
function of angular offset from the primary gamma-ray direction. The 
image templates shown here are evaluated at an impact distance (p ) of 
150 m and a first interaction depth (A) of 1 Xo. 

we note that the likelihood reconstruction can also be 
applied using templates generated with semi-analytic 
shower models l36l . 

Because the templates are produced from a simula¬ 
tion of the shower, the image model incorporates all ef¬ 
fects that influence the measured image intensity includ¬ 
ing atmospheric attenuation, geomagnetic field, tele¬ 
scope optics, and telescope detector response. The im¬ 
age model is a continuous distribution for the shower 
photons in the focal plane and the same template can 
therefore be used to compute the image intensity for 
cameras with arbitrary pixel geometry and field-of- 
view. For this study, we use image templates computed 
for the baseline telescope model with D = 12 m and 
A opt = 11.18 m 2 . The image intensity for other tele¬ 
scopes is calculated by rescaling the image intensity 
by the ratio of the telescope light collection area to the 
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baseline telescope model. 

The image intensity templates are stored on a six¬ 
dimensional grid: 

• log I0 (Energy) and Interaction Depth (log [(l /'’ and 
A) 

• Core Impact Distance and Position Angle (p and <p) 

• Projected Zenith and Azimuth Offset in Image 
Template Coordinates 

The shower templates are defined in a coordinate system 
rotated by <pi with respect to the shower axis such that 
the x-axis is aligned with the shower axis. In order to 
keep the memory footprint of the full six-dimensional 
template at a manageable level, the image templates are 
only defined over an angular region within one degree 
of the shower axis. 

The expected image intensity is computed from the 
image template sequence by a linear interpolation in the 
six-dimensional template space. The templates are also 
used to derive first derivatives of the image intensity as 
a function of the shower parameters which are used for 
calculation of the likelihood gradient. 

Figure [6] shows the image templates evaluated for 
gamma-ray showers of three different energies. The 
primary energy affects both the total intensity of the 
shower image as well as its shape. Higher energy show¬ 
ers propagate further into the atmosphere and result in 
shower images that are more extended along the shower 
axis. The core impact distance sets the geometrical 
perspective of the telescope and selects the Cherenkov 
light emission from particles with a specific range of 
angles with respect to the shower axis. Showers ob¬ 
served inside the Cherenkov light pool (p < 120 m) 
appear both brighter and narrower as the telescope ac¬ 
cepts Cherenkov light from higher energy particles that 
are closely aligned with the shower primary. More dis¬ 
tant showers are dimmer and increasingly offset from 
the primary origin. The interaction depth sets the start¬ 
ing point of the shower and primarily influences the dis¬ 
placement of the shower image along the shower axis. 

In the absence of the GF, the shower template is sym¬ 
metric with respect to the shower position angle. The 
GF breaks the axial symmetry as the Lorentz force pref¬ 
erentially perturbs the trajectory of the shower particles 
into the plane orthogonal to B . The GF effect is es¬ 
pecially pronounced for showers with small interaction 
depth for which the average propagation distance be¬ 
tween the first and second interactions is large. Figure 
[Tjillustrates the impact of the GF on the image template 
for three values of the core position angle ((f)): 0° (par¬ 
allel), 45°, and 90° (perpendicular). The shower width 


monotonically decreases as the shower position angle is 
increased from 0° to 90° reflecting the asymmetry in the 
shower development. For intermediate viewing angles 
(<p = 45°), the GF also causes a rotation of the image 
major axis relative to the shower axis. 


3.3. Gamma/Hadron Separation and Cut Optimization 

The final stage of the event analysis determines 
parameters that can be used for discrimination be¬ 
tween cosmic- and gamma-ray initiated air showers. 
The Cherenkov images produced by cosmic-ray show¬ 
ers can generally be distinguished from gamma-ray 
showers by their wider and more irregular appearance. 
Hadronic subshowers may also produce isolated clus¬ 
ters of Cherenkov light in the telescope image plane. A 
widely used set of parameters for background discrim¬ 
ination are the so-called mean scaled parameters de¬ 
fined already in the HEGRA collaboration If34l and ex¬ 
tensively used by current generation IACT experiments 
(see e.g., l37l ). The mean scaled parameters provide 
a measure of the deviation between the observed and 
expected telescope image moments for a gamma-ray 
shower. Using a set of simulated gamma-ray showers, 
lookup tables for the mean and standard deviation of the 
image moment parameters are produced as a function of 
the telescope image size and telescope impact distance 
(denoted here as p(S,p ) and cr p (S,p)). For a telescope 
image parameter p, we define the array-level parameter 
as 



Z 

i 


Pi - p(Si.pi) 

Wi -7T- 

o- p (Si,pi) 


( 21 ) 


where the sum is over all telescopes with reconstructed 
image parameters and vv, is a weighting factor. We use 
Wi = S i/cr p (S i,pi) assigning a larger weight to tele¬ 
scopes with brighter images and a smaller expected dis¬ 
persion in the image parameter. 

A second class of discriminant variables can be ob¬ 
tained by computing a goodness-of-fit between the data 
and the image template model evaluated at the best-fit 
shower parameters |[36l . When considering Gaussian- 
distributed data the natural goodness-of-fit parameter is 
the x~ statistic. For the purposes of background dis¬ 
crimination, it is not critical to have an exact model for 
the asymptotic distribution of the test statistic as long as 
it provides good separation power between signal and 
background. To quantify the agreement between the 
measured and expected pixel signals we define a ^ 2 -like 
parameter which we call the goodness-of-fit. 
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Figure 7: Image intensity templates as a function of angular offset from the primary gamma-ray direction for three values of the shower position 
angle: (f> = 0°, (f) = 45°, 0 = 90°. The templates shown are evaluated for a gamma-ray shower with an energy of 100 GeV, an impact distance of 
150 m, and an interaction depth of 0.3 Xq. The solid white line in each image shows the projection of the primary trajectory to the image plane. 
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■p. 




Pu{0) + Pb 


( 22 ) 


where Vi is a set of pixels in telescope i and N is the 
total number of pixels in the summation. We found that 
the best separating power was achieved by evaluating 


Equation 22 using the set of telescope pixels that sur¬ 
vive image cleaning, which we refer to as the image 
goodness-of-fit. 

To maximally exploit the rejection power drawn from 
the ensemble of event parameters we further make use 
of boosted decision trees (BDTs) generated with the 
TMVA package 081 . The use of machine learning tech¬ 
niques have been shown to provide significant improve¬ 
ment in overall background rejection power when ap¬ 
plied to IACT data {39]. We specifically use BDTs 
trained with the GradBoost algorithm with 200 trees, a 
maximum depth of 8, and a shrinkage parameter of 0.1. 
We train the decision tree analysis using the following 
six parameters: mean scaled width, mean scaled length, 
mean scaled displacement, array core distance, first in¬ 
teraction depth, and image goodness-of-fit. In order to 
avoid overtraining we use a training data set that con¬ 
stitutes 20% of the total gamma-ray and proton Monte 
Carlo samples. 


4. Results 


Using the simulation and analysis frameworks de¬ 
scribed in Sections [2] and [3] we have explored the in¬ 
fluence of the telescope design on the sensitivity and 
gamma-ray reconstruction performance of various array 
design concepts. Section 4.1 outlines the performance 
metrics used for comparison of the arrays. Section [4~2] 
defines a reference array alongside several benchmark 


configurations which are representative of realistic tele¬ 
scope and array configurations that will be chosen for 
CTA. In the subsequent sections we examine the influ¬ 
ence of each telescope parameter on the array perfor¬ 
mance. In Section 4.12 we study the performance of the 
benchmark arrays. 


4.1. Performance Metrics and Cut Optimization 

Our primary metric for the comparison of different ar¬ 
ray and telescope designs is the differential gamma-ray 
point-source sensitivity evaluated following the stan¬ 
dard procedure for CTA-related studies IT81 . The dif¬ 
ferential sensitivity is evaluated in a sequence of loga¬ 
rithmic bins of reconstructed energy with a width of 0.2 
dex (five bins per decade of energy). In each energy bin 
we calculate the expected number of signal and back¬ 
ground events within an energy-dependent aperture of 
radius 6. The number of signal events is estimated as¬ 
suming a gamma-ray point-source in the center of the 
FoV. The residual background rate is estimated by scal¬ 
ing the number of background events reconstructed in 
the inner 3° of the camera to the gamma-ray extrac¬ 
tion area. In each bin we require a 5<x excess above 
background and at least 10 signal events. The source 
significance is calculated using the method of m and 
a signal-free background region with a solid angle five 
times larger than the signal aperture. We further require 
a signal with a fractional amplitude above background 
of 5% in order to account for systematic errors in back¬ 
ground estimation. 

Sensitivity to spatially extended gamma-ray sources 
is calculated following the same procedure but with the 
gamma-ray signal spread out uniformly over a disk with 
angular diameter D. For a source with a given flux, the 
diffuse source sensitivity is always worse than the point- 
source sensitivity. In the case of a point-source, the 
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sensitivity depends on both background rejection effi¬ 
ciency and the PSF. The diffuse-source sensitivity, how¬ 
ever, depends primarily on the background rejection ef¬ 
ficiency and is nearly independent of the PSF when D is 
larger than the PSF. 

The quality of the gamma-ray reconstruction is esti¬ 
mated from the simulated gamma-ray PSF, shower core 
resolution, and energy resolution. The most important 
of these quantities is the gamma-ray PSF as it directly 
impacts the sensitivity to point-sources and the mor¬ 
phology of extended gamma-ray sources. We charac¬ 
terize the gamma-ray PSF by the radius that contains 
68% of the distribution of reconstruction errors (68% 
containment radius). 

When evaluating the performance of an array we ap¬ 
ply several selection criteria to reject both background 
events and gamma-ray events with poor reconstruction 
quality. Point-source cuts are composed of two energy- 
dependent selections on the gamma/hadron rejection pa¬ 
rameter and aperture radius, £{E) and 6(E), parame¬ 
terized as a cubic spline. The shape of these param- 
eterizations is independently optimized for each array 
and exposure time to maximize the differential point- 
source sensitivity versus energy. At high energies where 
the sensitivity of IACT arrays transitions from being 
background- to signal-limited, the optimal point-source 
sensitivity is obtained by increasing the gamma-ray effi¬ 
ciency and including events with poorer reconstruction 
quality and a higher background contamination level. 
Diffuse-source cuts are used when evaluating diffuse 
source sensitivity and comprise the same selection on 
the gamma/hadron parameter but with the aperture size 
fixed to the radius of the source (6(E) = D/2). 

Reconstruction cuts are used to define a homoge¬ 
neous sample of well-reconstructed showers with core 
locations within a predefined fiducial area of the array. 
Showers passing reconstruction cuts must have an im¬ 
pact distance from the array center that is less than 1.2 
times the distance from the center of the array to the 
nearest point along the array edge (as defined by the 
outer ring of telescopes). Showers with core locations 
near or within the array boundary (contained events) 
are sampled by a large number of telescopes that view 
the shower from multiple perspectives and allow for a 
more precise stereoscopic reconstruction of the shower 
trajectory. In arrays with mean telescope separations 
on par with the Cherenkov light pool size, contained 
events will also have one or more telescopes that sample 
the shower within its Cherenkov light pool, where the 
Cherenkov light from the highest energy shower parti¬ 
cles is visible. The light emitted from these particles 
provides a much better constraint on the shower trajec- 


Table 1: Geometrical characteristics and optical performance of the 
camera and optical systems of the DC-MST, SC-MST, and LST tele¬ 
scope designs chosen for the prod-2 MC design study ED The off- 
axis PSF performance is evaluated at a field angle equal to 3/4 of the 
distance to the edge of the FoV. For the LST and DC-MST designs 
which have cameras composed of hexagonal pixels, the given value 
of Dpj x is the width of a square pixel with the same solid angle as the 
hexagonal pixel used in that design. 



LST 

SC-MST 

DC-MST 

£pix [°] 

0.084 

0.067 

0.171 

A m [m 2 ] 

412 

50 

100 

A opt [m 2 ] 

52.5 

7.29 

13.65 

flpsf [°] (on-axis) 

0.03 

0.04 

0.04 

flpsf [°] (off-axis) 

0.12 

0.04 

0.08 

FoV [°] 

4.5 

8 

8 


tory than the light emitted by lower energy shower par¬ 
ticles. Events outside the array boundary (uncontained 
events) are sampled by a smaller number of telescopes 
for which the viewing angles are more closely aligned. 
This results in a less precise determination of the shower 
trajectory. 

Reconstruction cuts provide a measure of the gamma- 
ray reconstruction performance that can be evaluated in¬ 
dependently of the source strength and exposure time. 
Relative to point-source cuts, reconstruction cuts offer 
worse point-source sensitivity but a significantly better 
gamma-ray PSF at high energies (above 1 TeV). The 
improvement in the gamma-ray PSF can be attributed 
to the removal of uncontained events which are bright 
enough to trigger the array at high energies. This se¬ 
lection is very useful when studying strong sources to 
check morphology and spectral features while not rely¬ 
ing on the best signal-to-noise ratio. 

4.2. Benchmark Arrays 

The baseline CTA concept is an array of 50-100 tele¬ 
scopes distributed over an area of ~4-5 km 2 and com¬ 
posed of small-, medium-, and large-sized telescopes. 
Previous simulation studies have found that intermedi¬ 
ate layouts with a few (3-4) large-sized, on the order of 
20 medium-sized, and 50-60 small-sized telescopes of¬ 
fers the best tradeoff in performance over the full energy 
range of CTA 1(181 . Table □ shows the primary char¬ 
acteristics of the currently considered designs for the 
large- and medium-sized telescopes. The LST design 
is optimized for sensitivity at gamma-ray energies be¬ 
low 100 GeV and features a large effective mirror area 
which enables efficient triggering and reconstruction of 
low energy showers. The DC-MST and SC-MST are 
two alternative designs for the medium-sized telescope 
that would provide sensitivity in the core energy range 
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of CTA (100 GeV-10 TeV). The DC-MST is a single 
dish telescope that is similar in overall design to current 
generation IACTs with a camera composed of hexago¬ 
nal pixels with flat-to-flat spacing of 0.184°. [^] The SC- 
MST employs the dual-mirror Schwarzschild-Couder 
optical design and uses a smaller pixel size of 0.067° 
to achieve higher resolution imaging of the gamma-ray 
shower. Ray-tracing simulations of the SC-MST OS 
with realistic alignment tolerances for the mirrors and 
camera focal plane have demonstrated an optical PSF 
with a 68% containment radius between 0.02° and 0.04° 
ED- Although the array designs considered for previ¬ 
ous MC studies only incorporated DC-MSTs, the higher 
angular resolution SC-MST provides a potentially com¬ 
pelling option for the medium-sized CTA telescope. 

We consider several different benchmark array con¬ 
figurations shown in Table[2]to explore the performance 
of different array and telescope designs for CTA. M61 
is a reference array configuration with an effective light 
collection area that is intermediate between the DC- and 
SC-MST designs and an imaging performance similar 
to the SC-MST. M61SC is an array configuration with 
the same imaging performance as M61 but with a re¬ 
duced light collection area that is more comparable to 
the SC-MST design. M61DC and M25DC are chosen to 
represent a 61 and 25 telescope array respectively com¬ 
posed of telescopes of the DC-MST design. The latter 
configuration corresponds to the number of MSTs in the 
baseline CTA design 1181 . Arrays L5 and L61 are com¬ 
posed of telescopes with an optical effective area com¬ 
parable to the LST design and an imaging performance 
similar to the SC-MST. 

We show in Figure [8] the trigger effective area for ar¬ 
rays M61, M61SC, M61DC, and L61. The camera trig¬ 
ger threshold of each array is set using the telescope ef¬ 
fective light collection area and Equation [4] The sharp 
downturn in the effective area of the MST arrays around 
100 GeV can be attributed to the onset of the trigger 
energy threshold. Below the trigger threshold energy, 
the image amplitude of an average gamma-ray shower 
is insufficient to trigger telescopes within the Cherenkov 
light pool. At these energies only showers with large 
interaction depth can be effectively recorded, and the 
total effective area is primarily determined by the trig¬ 
ger efficiency for contained showers that impact within 
the array perimeter. At higher energies all of the arrays 
become fully efficient for triggering contained showers 
and differences in the effective area arise predominantly 
from the efficiency for detecting showers around the ar¬ 
ray perimeter. As the shower energy increases, the area 


3 Equivalent in solid angle to square pixels with D plx = 0.171°. 



Figure 8: Trigger effective area versus gamma-ray energy for arrays 
M61SC, M61. M61DC, and L61. The camera trigger thresholds (7’ t |,) 
for these arrays are 45, 60, 80, and 123 PE respectively. 

over which the arrays are fully efficient continues to in¬ 
crease and eventually extends well beyond the physical 
footprint of the array. Relative differences in the effec¬ 
tive area for telescopes with different A opt are signifi¬ 
cantly smaller at the highest energies as gains in the ef¬ 
fective area only come from showers around the array 
perimeter. 

4.3. Comparison with Other Telescope Simulations 

The simplifications in the detector response of the 
FAST simulation yield a much faster simulation code 
and enables the study of a broader parameter space 
compared to more detailed telescope simulations. Our 
simplified telescope model also enables us to employ 
a shower likelihood model which is nearly perfectly 
matched to the response characteristics of the tele¬ 
scopes. These simplifications allow us to explore the 
theoretical limit of the performance achievable by an 
1ACT array when all characteristics of the telescope op¬ 
tics and camera are accounted for in the event recon¬ 
struction. 

We have assessed the impact of the simplifications 
made in the FAST tool on the derived point-source 
sensitivity and gamma-ray PSF by comparing FAST 
against the well-tested sim_telarray package. We 
use both packages to simulate a 61 telescope array 
with the same geometry as our benchmark array lay¬ 
out with 120 m inter-telescope separation. For the 
sim_telarray simulation we use the prod-2 SCT 
model lf33l with a trigger pixel threshold of 3.1 PE. For 
the FAST simulations, we use a telescope model with 
the same performance characteristics as the prod-2 SCT 
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Table 2: Number of telescopes and telescope model parameters of the benchmark array configurations used for this study. All arrays are composed 
of telescopes arranged on a uniform grid with constant inter-telescope spacing of 120 m as shown in Figure^ 


Name 

N te l 

A op t [m 2 ] 

Dpi X n 

^psf [°] 

T th [PE] 

R nsb [MHz] 

M61 

61 

11.18 

0.06 

0.02 

60 

14.7 

M61SC 

61 

8.38 

0.06 

0.02 

45 

11.1 

M61DC 

61 

14.91 

0.16 

0.08 

80 

139.5 

M25DC 

25 

14.91 

0.16 

0.08 

80 

139.5 

F5 

5 

47.15 

0.06 

0.02 

123 

61.9 

F61 

61 

47.15 

0.06 

0.02 

123 

61.9 




10 2 10 3 10 4 
Energy [GeV] 


Figure 9: Performance of a 61 telescope array simulated with FAST (black squares) and sim.telarray (blue circles). Left: Differential point- 
source sensitivity for a 50 h observation time. Shown as the solid gray line is the differential sensitivity of Array I from cm evaluated with the 
most sensitive analysis at each energy from the four alternative analyses presented in that work (MPIK, IFAE, SAM, and Paris-MVA). Right: 68% 
containment radius of the gamma-ray PSF after point-source cuts. 


model shown in Table [I] and a camera trigger threshold 
of 42 PE. Relative to the telescope model used for the 
M61SC benchmark array, the prod-2 SCT model has a 
larger pixel size and optical point-spread function and 
a slightly smaller light collection area. For gamma-ray 
showers near the trigger threshold (E ~ 100 GeV), the 
sim_telarray telescope model has a slightly lower ef¬ 
fective camera threshold as compared to the telescope 
in our simplified simulations. The choice of a higher 
threshold for our simulations was made to be conser¬ 
vative and limit possible overestimations in sensitivity 
close to the threshold. 

Fig. [9] shows the comparison of the array perfor¬ 
mance obtained when simulating the same array with 
sim_telarray and FAST. We include in the same fig¬ 
ure the point-source sensitivity of Array I from m 
which was simulated using sim_telarray but with a 
different array and telescope setup. At energies above 
75 GeV the point-source sensitivity obtained with the 
FAST simulations is 20% better than the sim_telarray 
simulations. The gamma-ray PSF and gamma-ray re¬ 


construction efficiency is similar over the same energy 
range indicating that the improvement in point-source 
sensitivity can be attributed to an enhanced gamma- 
hadron separation in the FAST simulations. Below 
50 GeV the sim_telarray simulations yield a 40% 
better sensitivity due to the slightly lower telescope 
trigger threshold. Although we observe measurable 
differences in the array performance when comparing 
our simulations with sim_telarray, the differences in 
point-source sensitivity are much smaller than the dif¬ 
ferences between individual analysis packages that use 
the same sim_telarray simulations as input (see e.g. 
the comparison of alternative analyses in Ifl8l ). Further¬ 
more the conclusions drawn in this work about the rel¬ 
ative performance of different array and telescope de¬ 
signs is most likely not affected by these differences. 
It is thus easy to scale our results and readily compare 
them to other sim_telarray results. 
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Figure 10: Reconstruction performance and gamma-ray point-source sensitivity of array M61 obtained with different event reconstruction and 
analysis algorithms: likelihood (black, solid), likelihood without goodness-of-fit (blue, dashed), and moment (red, dot-dashed). Left: Differential 
point-source sensitivity for a 50 h observation time. Right: 68% containment radius of the gamma-ray PSF after point-source cuts. 


4.4. Performance of the Likelihood Reconstruction 


Relative to moment-based reconstruction techniques, 
likelihood-based reconstruction algorithms have been 
shown to provide better gamma-ray angular resolution 
as well as improved separation between gamma-ray 
and cosmic-ray induced showers ED ESI ED. We as¬ 
sess the relative improvement from the likelihood ap¬ 
proach by comparing its performance with an analysis 
that uses only the geometric trajectory reconstruction 
and moment-based image parameterization described in 
Section 13.21 which we refer to here as the moment re¬ 
construction. Because the moment reconstruction is 
more sensitive to the presence of noise fluctuations in 
the image, we use a slightly higher cleaning thresh¬ 
old (s/cr = 9) than the threshold used for the likeli¬ 
hood analysis. We use a BDT background discriminant 
trained with the same settings described in Section 3.3 
but excluding parameters derived from the likelihood 
analysis. 

Figure [10] shows the comparison of the point-source 
sensitivity obtained with the moment analysis, the like¬ 
lihood analysis, and a likelihood analysis in which the 
goodness-of-fit (GOF) parameter is excluded from the 
training of the decision tree. With the likelihood-based 
analysis, we find a factor of two improvement in point- 
source sensitivity and a 30-40% improvement in the 
gamma-ray PSF over the full energy range. As seen 
from the comparison between the likelihood analyses 
performed with and without the GOF parameter, the im¬ 
provement in the point-source sensitivity is attributable 
to gains in both the gamma-ray angular resolution and 
the background rejection power. The addition of the 


GOF parameter provides an additional 30% improve¬ 
ment in sensitivity. 

We also observe that the energy threshold for the like¬ 
lihood analysis is considerably lower (Ifu - 50 GeV) 
relative to the moment reconstruction (E^ - 100 GeV). 
The improved performance of the likelihood analysis at 
low energies can be attributed to both the higher im¬ 
age reconstruction efficiency and the smaller bias of the 
likelihood energy estimator. Because the likelihood re¬ 
construction is insensitive to the inclusion of pixels with 
small signals, the cleaning threshold can be optimized to 
maximize the reconstruction efficiency for low-energy 
showers without impacting the performance at higher 
energies. 

4.5. Influence of the Geomagnetic Field 

The deflection of charged particles in the EM shower 
by the geomagnetic field (GF) can significantly distort 
the shapes of gamma-ray images recorded by IACTs 
El 03 04). The strength and orientation of the GF is 
thus an important consideration for the selection of can¬ 
didate sites for an IACT observatory. Its influence can 
be as large or larger than the site elevation 1145)1 . The 
magnitude of the induced deflection is proportional to 
the perpendicular component of the GF ( B ± ) and there¬ 
fore the strength of the GF effect depends on both the 
magnitude of the GF vector as well as its orientation 
relative to the shower trajectory. Due to the asymmetry 
in the shower shape induced by the GF, the distortion 
visible to a telescope also depends on the orientation of 
the shower impact point relative to the telescope posi¬ 
tion. Telescopes with shower position angles close to 0° 
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Figure 11: Performance of array M61 simulated with the equatorial GF (B ± = 20.7 yuT; red diamonds and solid line), a GF configuration with 
a reduced perpendicular component (B± = 10.35 f/T; blue circles and solid line), and no GF (black squares and solid line). Left: Differential 
point-source sensitivity for a 50 h observation time. Right: Gamma-ray angular resolution (68% containment radius) after reconstruction cuts. 
Dashed curves show the same comparison for gamma-ray showers with an interaction depth (i) greater than 1.0 Xq. 


or 180° see a larger GF effect as the GF-induced elon¬ 
gation in the shower occurs primarily in the plane per¬ 
pendicular to the telescope pointing. 

To obtain a realistic assessment of the GF effect for 
any given observatory site would require simulations 
with many telescope orientations as they occur for re¬ 
alistic observation profiles of gamma-ray sources. We 
did not carry out such simulations and instead focused 
on the effect of the GF for a few representative values 
of B ± . Our baseline site configuration has (B X ,B Z ) = 
(27.5 yuT, -15.0 yuT) with B ± = 20.7 /ZT when observ¬ 
ing a shower with Zn = 20° and Az = 0°. To test the 
influence of the GF strength we performed simulations 
of array M61 for two additional site models: a site with 
(B x , B -) = (19.84 /rT, -24.24 //T) that has a perpendicu¬ 
lar GF component that is half as large as for our baseline 
site (B ± = 10.35 //T) and a site with no geomagnetic 
field. By limiting ourselves to these few cases we can 
only give a general guidance for effects of the magnetic 
field on any observable gamma-ray source. Depending 
on the specific source observation profile, the effect of 
the GF for an individual source might be different. 

The configurations we tested have a range of 
field strengths that are comparable to the southern 
Hemisphere sites considered for CTA. The Namib¬ 
ian H.E.S.S. site and the Argentinian Leoncito sites 
have (B X ,B Z ) = (12.1 yuT,-25.5 /tT) and (B x , B z ) = 
(20.1 //T,-12.2 /ZT), respectively [[431 . Because the 
strength and orientation of the GF is generally a slowly 
varying function of the site latitude and longitude these 
two sites provide a good representation of the expected 


GF effect for sites in Africa and South America. When 
observing a shower at Zn = 20° and Az = 0° the Namib¬ 
ian and Argentinian sites have perpendicular compo¬ 
nents of 2.7 /ZT and 14.7 /ZT. However a more realistic 
measure of the expected GF effect is the average perpen¬ 
dicular component over the range of azimuth angles that 
a gamma-ray source is observed. The Namibian and Ar¬ 
gentinian sites have an average GF strength at Zn = 20° 
of 13.4 /ZT and 19.7 /ZT, respectively. 


The comparison of the array performance for the 
three GF configurations is presented in Figure [IT] We 
find that the effect of the GF strength is strongest at 
100 GeV where the point-source sensitivity is reduced 
by 50% when increasing B_ from 0 //I’ to 20.7 /IT. We 
also observe that the effect of the GF scales linearly with 
B ± such that the site configuration with B ± = 10.35 /rT 
suffers approximately half of the reduction in sensitiv¬ 
ity relative to our baseline site configuration. Below 
energies of 100 GeV, the effect of the GF is lessened 
because only gamma rays that convert deep in the at¬ 
mosphere can be efficiently reconstructed. The lower 
the particle interacts in the atmosphere the less it is af¬ 
fected by the GF. At higher gamma-ray energies the im¬ 
pact of the GF is lessened due to both the higher energy 
of the secondary particles and the larger path length in 
the atmosphere. As shown in the right panel of Fig¬ 
ure 11 the GF worsens the point-source sensitivity pri¬ 


marily by degrading the gamma-ray PSF. For showers 
with interaction depths larger than 1 Xo, differences in 
the gamma-ray PSF between the different GF configu¬ 
rations are found to be less than 20% illustrating that 
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the influence of the GF increases with decreasing inter¬ 
action depth. 


4.6. Imaging Performance 

The telescope design has a large impact on the result¬ 
ing gamma-ray PSF obtained with the complete array. 
The optical design of the individual telescopes defines 
their achievable optical PSF and the camera design de¬ 
termines how efficiently the optical PSF can be trans¬ 
lated into an improved gamma-ray PSF. For a given op¬ 
tical PSF, the gamma-ray PSF can be improved by re¬ 
ducing the camera pixel size. In the limit that the pixel 
size is much smaller than the PSF, the improvement of 
the gamma-ray PSF saturates and a further reduction in 
pixel size does not provide any measurable advantage 
but increases the cost of the camera. Thus the optimal 
tradeoff between performance and cost is one in which 
the pixel size is appropriately matched to the quality of 
the optical PSF. Current generation IACTs have cameras 
using pixel sizes from 0.1° to 0.16° and an optical PSF 
at the center of the FoV which is considerably smaller 
than the pixel size. Here we explore a new parameter 
space for the IACT imaging resolution by examining the 
performance of camera designs with pixel sizes between 
0.04° and 0.1°. Such designs begin to resolve the core 
of the Cherenkov shower which has an intrinsic angular 
size of ~0.01°. 

The left panel of Fig. 12 shows the gamma-ray PSF 
versus pixel size for arrays with different optical PSFs. 
For an optical PSF of 0.08° the gamma-ray PSF shows 
only a modest improvement of ~ 10% when reducing 
the pixel size from 0.2° to 0.04°. An optical PSF be¬ 
tween 0.02° and 0.04° is found to be critical to real¬ 
ize the full improvement in gamma-ray PSF that can be 
achieved by reducing the camera pixel size below 0.12°. 
The improvement of the gamma-ray PSF at different en¬ 


ergies when reducing the pixel size is shown in Fig. 12 


The gamma-ray PSF is significantly better at all ener¬ 
gies when reducing the pixel size. There is a slight mod¬ 
ulation seen in the improvement versus energy more 
pronounced for larger pixel sizes. The smaller pixel size 
performs best at low and high energies (E < 100 GeV 
and E > 2.5 TeV) while the improvement is less pro¬ 
nounced in the intermediate energy range. An improve¬ 
ment of the gamma-ray PSF of about 20% in the full 
energy range by reducing the pixel diameter from 0.12° 
to 0.06° demonstrates a realistic difference between cur¬ 
rently considered optical telescope designs for CTA. 

The effect of the pixel size on the differential point- 
source sensitivity is shown in FigJT3] The pixel size has 
the strongest impact at low energies (< 100 GeV) where 
a factor of two relative improvement is observed when 


the pixel size is reduced from 0.16° to 0.06°. At higher 
energies the smaller pixel size results in a smaller but 
still measurable improvement in point-source sensitiv¬ 
ity of 30-40%. Above 3 TeV differences between ad¬ 
jacent pixel sizes become indistinguishable due to the 
limited background statistics that make evaluation of 
small sensitivity differences very difficult. The gamma- 
ray PSF is clearly improved over the complete energy 
range by about 50% as the pixel size is reduced from 
0.2° to 0.06°. The observed improvement in sensitiv¬ 
ity demonstrates that the intrinsic shower features that 
can be used for background suppression and direction 
reconstruction are still smaller than the pixel sizes of 
currently operating Cherenkov telescopes. 


4.7. Light Collection Area 


The telescope light collection area determines the 
signal-to-noise ratio (SNR) of the shower images and 
the efficiency with which these images can be recorded 
by the telescope trigger. Therefore we expect that a 
larger A opt increases the trigger efficiency and provides 
better defined images and hence improves performance 
of the array. The role of the A opt parameter is partic¬ 
ularly relevant for the performance of the array at low 
energies where the smaller light yield per image makes 
reconstruction and analysis of the gamma-ray showers 
more challenging. 

The assumed design, size, and cost of the proposed 
telescopes yields distinct A opt values. We studied the ef¬ 
fect of the A opt on the gamma-ray PSF and point-source 
sensitivity of the array by examining the performance 
of telescope models with A opt between 2 m 2 and 50 m 2 . 
These models span the range of light collection areas 
between SST-like and LST-like telescope designs. The 
SST, MST, and LST telescope designs have A opt of ap¬ 
proximately 1-2 m 2 , 5-10 m 2 , and ~50 m 2 respectively 

EDESII. 


Figure 14 shows the comparison of the gamma-ray 
PSF and point-source sensitivity for telescopes with A opt 
between 1.98 m 2 (SST-like) and 47.15 m 2 (LST-like). 
A opt has only a minor effect on the gamma-ray PSF in 
most of the energy range investigated here. In the mid¬ 
dle energy range between 100 GeV and 1 TeV we find 
an improvement of 5-10 % when increasing the tele¬ 
scope light collection area from 11.18 m 2 to 47.15 m 2 . 
The almost insignificant improvement around 100 GeV 
is caused by a selection effect of the reconstructed 
gamma-ray events. At these low energies, telescopes 
with smaller A opt can only trigger on the brightest show¬ 
ers that convert deep in the atmosphere. As discussed in 
Section 4.5 the larger interaction depth of these showers 
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Figure 12: Left: 68% containment angle of the gamma-ray PSF at 317 GeV versus camera pixel size for telescope models with different optical 
PSFs (/? ps f): 0.02° (black squares), 0.04° (blue circles), 0.08° (red diamonds). The gamma-ray PSF is evaluated after applying reconstruction cuts. 
The baseline configuration for all simulations is array M61. Right: 68% containment angle of the gamma-ray PSF versus gamma-ray energy for 
array M61 with R ps f = 0.02° simulated with different telescope pixel sizes: 0.06° (black squares), 0.12° (blue circles), 0.20° (red diamonds). 


lessens the impact of the GF and results in a more ac¬ 
curate reconstruction of the direction. Larger telescopes 
can efficiently trigger on showers with both large and 
small interaction depths which results in a larger effec¬ 
tive area but a worsening of the overall gamma-ray PSF. 
This effect reverses at the very lowest energies (30- 
50 GeV) where the reduced SNR images recorded by 
telescopes with small A opt dominates the reconstruction 
quality. 

The light collection area has a measurable impact 
on the point-source sensitivity only at energies below 
300 GeV with telescopes with larger light collection 
area yielding better sensitivity. The increase in sensi¬ 
tivity is most significant below 100 GeV and is a re¬ 
sult of the reduction in the telescope trigger threshold 
and resulting increase in the gamma-ray effective area. 
The larger light collection area also yields better SNR in 
the shower images improving the reconstruction of low 
energy events. At higher energies the impact of light 
collection area is significantly reduced as the array be¬ 
comes fully efficient for triggering and reconstructing 
events that impact within the array boundary. Improving 
the image SNR provides little improvement at these en¬ 
ergies because the reconstruction is predominantly lim¬ 
ited by intrinsic shower fluctuations. Remarkably the 
improvement in point-source sensitivity is almost neg¬ 
ligible between telescopes with 26.51 m 2 and 47.15 m 2 
over the whole energy range. 

The observed improvements in array performance 
above the trigger threshold are small when consider¬ 
ing that light collection area is the dominant parameter 


influencing the telescope cost. Given the small differ¬ 
ences in reconstruction performance, the primary moti¬ 
vation for choosing a telescope design with larger light 
collection area is to reduce the array energy threshold. 
However for an array of fixed cost increasing the light 
collection area also entails a reduction in the number of 
telescopes. For gamma-ray energies between 100 GeV 
and 1 TeV, a telescope with A opt of 5-10 m 2 (MST-like) 
clearly provides the best performance to cost ratio. Ar¬ 
ray designs that include a small number of telescopes 
with larger light collection area can lower the energy 
threshold while keeping the cost of the total array within 
reasonable limits. Performance of arrays with differ¬ 
ent numbers of telescopes are studied further in Section 

ELD 

4.8. Inter-telescope Separation 

The inter-telescope separation determines both the 
physical area of the array footprint as well as the average 
number of telescopes that will participate in the recon¬ 
struction of individual showers. Smaller telescope sep¬ 
arations improve reconstruction quality for contained 
showers at the cost of lowering the total effective area 
of the array. Larger telescope separations are generally 
preferred when optimizing for sensitivity at higher en¬ 
ergies since the point-source sensitivity of IACT arrays 
at moderate exposures (10-50 hours) is signal limited 
above 1-3 TeV. Another important consideration when 
optimizing the telescope separation is the number of 
telescopes within the Cherenkov light pool. Telescopes 
within the Cherenkov light pool sample light emitted by 
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Figure 13: Performance of array M61 simulated with pixel sizes from 0.04° to 0.20°. Left: Differential point-source sensitivity for a 50 h 
observation time. Right: 68% containment angle of the gamma-ray PSF evaluated after point-source cuts. 
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Figure 14: Performance of array M61 simulated with different values of A op t: 1.98 m 2 (black squares), 4.71 m 2 (blue circles), 11.18 m 2 (red 
diamonds), 26.51 m 2 (magenta triangles), and 47.15 m 2 (cyan triangles). Left: Differential point-source sensitivity for a 50 h observation time. 
Right: 68% containment angle of the gamma-ray PSF after reconstruction cuts. 


higher energy particles in the shower core and provide 
a more accurate determination of the shower trajectory. 
Telescope separations that are comparable to the size of 
the Cherenkov light pool (100-150 m) ensure that mul¬ 
tiple telescopes will sample each shower within its light 
pool. Finally smaller separations may potentially im¬ 
prove background rejection by increasing the efficiency 
for detecting Cherenkov light from hadronic subshow¬ 
ers produced in cosmic-ray background events. 


The impact of the telescope separation on the gamma- 
ray PSF is illustrated in the top panel of Fig. 15 which 
shows a comparison of arrays with separations between 
60 m and 200 m. In this comparison we consider only 
showers passing reconstruction cuts with core positions 


near or within the array boundary. These cuts select 
events with the best PSF and reduce the differences in 
performance caused by the finite array size. The re¬ 
duction of the telescope grid spacing from 120 m to 
60 m results in a 20% improvement of the gamma-ray 
PSF between 30 GeV and 10 TeV. However this rather 
small improvement would require a quadrupling in the 
number of telescopes to cover a similar area. Thus the 
improvement of the gamma-ray PSF from reducing the 
telescope spacing has to be compared to the reduction of 
effective detector area when fixing the number of avail¬ 
able telescopes. 

The lower left and right panels of Fig. [l5]show the 
gamma-ray PSF and point-source sensitivity for the set 
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Figure 15: Performance of array M61 simulated with different inter-telescope separations: 60 m (black squares), 80 m (blue circles), 120 m (red 
diamonds), 160 m (magenta triangles) and 200 m (cyan triangles). Top Left: 68% containment angle of the gamma-ray PSF after reconstruction 
cuts. Top Right: Gamma-ray effective area after point-source cuts. Bottom Left: Differential point-source sensitivity for a 50 h observation time. 
Bottom Right: 68% containment angle of the gamma-ray PSF after point-source cuts. 


of telescope separations evaluated with a selection op¬ 
timized for point-source sensitivity. The increase of 
effective area with larger telescope spacing generally 
outweighs the reduction of sensitivity due to a worsen¬ 
ing of the gamma-ray PSF. The point-source sensitivity 
improves with increasing telescope spacing at energies 
above 100 GeV with the best sensitivity achieved with 
a telescope spacing of 160-200 m. When increasing the 
telescope spacing to 200 m a noticeable worsening of 
the sensitivity below 300 GeV is seen because the num¬ 
ber of individual telescopes triggering on each event is 
reduced and hence the information available for direc¬ 
tion and particle type reconstruction. 


When evaluated with point-source cuts as shown in 
the bottom right panel of Fig. 15 the gamma-ray PSF 
above 300 GeV becomes worse as the telescope separa¬ 
tion is decreased. Although a smaller separation gives 


a better reconstruction for contained events, the smaller 
array footprint results in a larger fraction of uncontained 
events which tend to dominate the PSF at high energies. 
This emphasizes that for most applications where the 
maximum sensitivity of the array is required the PSF 
has a quite different behavior compared to the theoret¬ 
ically possible behavior. A wider spacing of the MSTs 
will provide a much better performance for most sci¬ 
ence cases compared to a narrow spacing that would be 
only beneficial for the very few cases where the gamma- 
ray PSF is much more important than sensitivity. Thus 
the best spacing for the MSTs for all purposes is about 
160 m. 

4.9. Trigger Threshold 

The telescope trigger threshold is an important quan¬ 
tity to determine the accessible energy range by any 
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Figure 16: Performance of array M61 simulated with different camera trigger thresholds: 34 PE (black squares), 45 PE (blue circles), 60 PE (red 
diamonds), 80 PE (magenta triangles). Left: Differential point-source sensitivity for a 50 h observation time. Right: Gamma-ray effective area 
after point-source cuts. 


telescope array. The impact of the individual tele¬ 
scope trigger threshold is studied on the differential 
point-source sensitivity of the M61 baseline array (see 
Fig. 16 1 . As expected for an MST-like array with A opt ^ 
10 rrr :he trigger threshold has little effect on the sen¬ 
sitivity at energies above 100 GeV. At higher energies 
the telescope trigger becomes fully efficient for show¬ 
ers impacting within the array and reducing the trigger 
threshold only increases the efficiency for showers on 
the array periphery. Because these distant showers are 
generally not well reconstructed they do not contribute 
to the array sensitivity. 


Reducing the telescope trigger threshold of Ar¬ 
ray M61 is found to significantly improve the point- 
source sensitivity below 100 GeV. A reduction of the 
trigger threshold from 80 PE to 34 PE results in a sig¬ 
nificant improvement at energies below 100 GeV and 
reaches up to an order of magnitude at 30 GeV. How¬ 
ever, in a realistic telescope design the accidental trigger 
rate can not be arbitrarily high due to the limitations on 
the readout rate that can be sustained by the telescope 
data acquisition. The 60 PE effective trigger threshold 
chosen for Array M61 is a realistic target for a trig¬ 
ger implementation that follows the same design used 
by current generation lACTs. Lower trigger thresholds 
may be achievable by employing more sophisticated de¬ 
signs for the camera- and array-level triggers such as re¬ 
quiring additional trigger topologies for individual tele¬ 
scopes or higher multiplicities for the array trigger. If 
further improvements in the performance of the trigger 
can be realized then the presented sensitivities at low en¬ 
ergies could be further improved. Furthermore, it is evi¬ 


dent that the likelihood reconstruction is very efficient at 
low energies and that any reduction in trigger threshold 
is directly translated into an improvement in sensitivity. 
The same statement is not necessarily true for the mo¬ 
ment reconstruction that usually has a higher analysis 
threshold compared to the likelihood reconstruction as 


shown in Fig. 10 


4.10. NSB Rate 

Night-sky background (NSB) is caused by the pres¬ 
ence of light sources such as stars, the Moon, and arti¬ 
ficial light pollution and represents an irreducible back¬ 
ground for the reconstruction and analysis of gamma- 
ray air showers. Because the Cherenkov photons de¬ 
tected in a single pixel have an intrinsic arrival time dis¬ 
persion of 3-6 ns, lACTs can significantly reduce the 
NSB by integrating the Cherenkov signal in a narrow 
time window (typically with AT ~ 10 ns). The inte¬ 
grated NSB level thus depends on both the NSB rate as 
well as the size of the window used for signal integra¬ 
tion. The need for a small integration window motivates 
camera designs with high bandwidth readout electronics 
which would allow the integration window to be made 
as small as possible. The impact of the NSB rate on the 
sensitivity of the array is also important when consid¬ 
ering possible observatory sites and performing obser¬ 
vations during moonlight. Moonlight observations can 
considerably increase the duty cycle of the observatory 
although the exact amount of observation time gained 
depends on the NSB rate that the individual telescope 
can handle. 

We studied the impact of NSB on the performance 
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Figure 17: Performance of arrays M61DC (red) and M61SC (blue) simulated with a baseline NSB flux of 365 MHz deg -2 m -2 (circles and solid 
lines) and an NSB flux that is 3 (dashed) and 6 (dash-dotted) times higher than the baseline value. Left: Differential point-source sensitivity for a 
50 h observation time. Right: Gamma-ray effective area after point-source cuts. 


of the array by performing simulations with three NSB 
flux levels: a baseline flux level with an integral flux 
of 365 MHz deg~ 2 m~ 2 and NSB fluxes that are 3 and 
6 times higher than the baseline flux. As described in 
Section |2.3| the baseline flux level corresponds to the 
expected night-sky intensity for a dark, extragalactic 
field. The higher NSB fluxes are representative of either 
a higher NSB rate due to operation under high night-sky 
brightness (moonlight) or a longer effective integration 
window. A higher NSB rate also increases the rate of 
accidental triggers and would require a higher trigger 
threshold in order to maintain the accidental trigger rate 
at a constant level. For this study we kept the trigger 
threshold fixed at its nominal value and only examine 
the impact of the NSB on the pixel SNR. 

Figure [IT] shows the comparison of the point-source 
sensitivity and gamma-ray effective area of arrays 
M61SC and M61DC simulated at the three NSB levels. 
The NSB level only appreciably affects the sensitivity 
below 300 GeV where the SNR of the shower image 
is lowest. Most of the reduction in sensitivity is a re¬ 
sult of the lower reconstruction efficiency as low SNR 
images are removed at the cleaning stage of the anal¬ 
ysis. Remarkably the reduction in sensitivity is much 
more pronounced in the case of larger pixels (DC-like 
telescope). In case of the SC telescope design, opera¬ 
tion at a six times higher NSB rate would only degrade 
the sensitivity below about 100 GeV and only up to a 
factor of two. The DC-like design would also suffer sig¬ 
nificant sensitivity loss only below about 100 GeV but 
to a much greater degree. Here it should be noted that 
the sensitivity advantage of the DC telescopes below 50 


GeV under low NSB is lost in case of three times in¬ 
creased NSB and that the SC design is better for six 
times higher NSB at all energies. 

4.11. Number of Telescopes in the Array 

One of the most important parameters concerning the 
sensitivity of an IACT array is the number of telescopes. 
A larger number of telescopes increases both the total 
effective area for triggering and reconstructing gamma- 
ray showers but also increases the average number of 
telescopes that participate in the reconstruction of each 
shower. Increasing the number of telescopes leads to 
better point-source sensitivity and an improved gamma- 
ray PSF. 

Figure [T8] compares the performance of arrays with 
between 5 and 61 telescopes. We investigate the scal¬ 
ing relation of the improvement in sensitivity with in¬ 
creasing number of telescope. In the limit of an infi¬ 
nite array the point-source sensitivity should scale with 
the number of telescopes as N^ 2 . However we observe 
an increase of sensitivity that is slightly better than the 
N^ 2 at all energies. This emphasizes that in the case 
of small telescope arrays increasing the number of tele¬ 
scopes yields larger improvements as compared to the 
case of extending large arrays. Adding 36 telescopes to 
a 25 telescope array improves the sensitivity by a factor 
of — 1.7-1.8. 

In contrast to the point-source sensitivity, the gamma- 
ray PSF improves non-uniformly over energy with in¬ 
creasing telescope number. The best improvement is 
seen at larger energies while at E < 300 GeV the im¬ 
provement is only clearly visible between 5 and 13 tele- 
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Figure 18: Performance of array layouts with telescope number (N te j) from 5 to 61. All arrays are simulated with a 120 m inter-telescope separation 
and the same telescope model as Array M61. Left: Differential point-source sensitivity for a 50 h observation time. Right: 68% containment angle 
of the gamma-ray PSF after applying point-source cuts. 


scopes. At high energies the curves in Fig. [18] show 
a clearer separation demonstrating that more telescopes 
help to better localize the showers above 1 TeV. The 
energy dependency has its origin in the fact that only 
high energy showers produce enough light to trigger dis¬ 
tant telescopes. Thus larger arrays with more telescopes 
benefit at high energies because the average number of 
telescopes participating in the shower reconstruction is 
increased. In the case of lower energy showers, the 
number of telescopes contributing to the shower anal¬ 
ysis is limited by the telescope spacing and not the ab¬ 
solute number of telescopes in the array. Increasing the 
footprint of the array also increases the parallax between 
telescopes observing an uncontained shower. The larger 
parallax yields a better shower direction reconstruction 
and further improves the reconstruction performance at 
high energies. 

4.12. Comparison of Array Designs for CTA 

After studying the effect of individual telescope pa¬ 
rameters on the point-source sensitivity and gamma- 
ray PSF, we now compare realistic telescope designs 
against each other to find a suitable array design for 
CTA. To achieve a comprehensive comparison we in¬ 
vestigate all the benchmark arrays defined in Table [2] 
and give a quantitative comparison between the differ¬ 
ent telescope layouts. Among the benchmark arrays are 
also two more theoretically interesting cases. Array L61 
is representative of the theoretical limit for an IACT ar¬ 
ray if the budget is not limited and only the number of 
telescopes is fixed. In a similar fashion. Array L5 is 
included to study the contribution of an LST subarray 


with 3-5 telescopes such as currently considered for the 
baseline configuration of CTA. 

Fig. [19] shows that Array M61SC is more sensi¬ 
tive than Array M61DC at all energies above 50 GeV, 
where the increase in sensitivity is about 30%. In ad¬ 
dition to the improvement in point-source sensitivity, 
the M61SC array also has a better gamma-ray PSF at 
all energies. The smaller gamma-ray PSF would help 
to determine the morphology of extended sources and 
help to separate point sources. These additional impor¬ 
tant effects are difficult to assess quantitatively because 
they heavily rely on the source population and proper¬ 
ties in the sky. The diffuse source is simulated as an 
uniformly extended disk with a radius of 0.5°. The 
diffuse-source sensitivity does not show any improve¬ 
ment of the M61SC array over the M61DC array be¬ 
cause the gamma-ray PSF does not help to reduce the 
background but still the M61SC would enable for a non- 
uniform source to asses the morphology better than Ar¬ 
ray M61DC. The diffuse source sensitivity emphasizes 
that the sensitivity gain of the SC array compared to the 
DC array comes almost entirely from the PSF improve¬ 
ment while the improvement in the background rejec¬ 
tion power is marginal. 

Array M25DC is representative of the MST subset 
of the CTA array design as it was planned without a 
US contribution. Comparing the Array M61SC and 
Array M61DC to the M25DC baseline configuration, 
it is obvious that adding MST telescopes will improve 
the sensitivity of CTA in the key energy range between 
100 GeV and about 1 TeV by about a factor two regard¬ 
less of their design. This is expected from the fact that 
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Figure 19: Performance of benchmark arrays: M61SC, M61DC, M25DC, L5, and L61. Top Left: 68% containment angle of the gamma-ray 
PSF after applying point-source cuts. Top Right: Gamma-ray effective area after point-source cuts. Middle Left: Differential rate of the total 
cosmic-ray background (protons and electrons) after point-source cuts. Middle Right: Differential rate of protons after point-source cuts. Bottom 
Left: Differential point-source sensitivity for a 50 h observation time. Shown as the solid gray line is the differential sensitivity of Array I from 
03 evaluated with the most sensitive analysis at each energy from the four alternative analyses presented in that work (MPIK, IFAE, SAM, and 
Paris-MVA). Bottom Right: Differential diffuse-source sensitivity (D = 0.5°) for a 50 h observation time. 
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the sensitivity is improved by the addition of telescopes, 


as shown in Fig. 18 


We also compared the point-source sensitivity of our 
benchmarks arrays with Array I from tm For Array I 
we use a differential sensitivity curve that is constructed 
by taking the best sensitivity in each energy bin from the 
four alternative analyses presented in that work (MPIK, 
IFAE, SAM, and Paris-MVA). Although the simulations 
in this paper were performed with different telescope 
models and a different detector simulation package, this 
array is representative of the expected performance of 
the baseline CTA concept. In the central energy range 
from 100 GeV to 3 TeV, Arrays M61DC and M61SC 
provide a factor of 3-4 improvement in point-source 
sensitivity relative to Array I. This improvement can 
be primarily attributed to the increase in the number of 
MSTs from 18 to 61. Array I performs better at ener¬ 
gies below 50 GeV and above 3 TeV as compared to 
Array M25DC and even Arrays M61DC and M61SC. 
This improvement can be attributed to the inclusion of 
56 SSTs and 3 LSTs in Array I. Array L5 was simu¬ 
lated with five LSTs very similar to the ones included 
in Array I, and the sensitivity curve obtained for L5 
matches very well the sensitivity of Array I at low ener¬ 
gies, demonstrating that the advantage of Array I at low 
energies does in fact come from the LSTs. 

Linally Array L61 yielded only an improvement be¬ 
low 100 GeV, making such an array impractical based 
on the large cost differential between a single MST and 
LST. However the performance of this array shows what 
is theoretically achievable in the case of no budget con¬ 
straints. Array M61SC provides comparable sensitivity 
to Array L61 at all energies above 100 GeV and thus is 
very close to the performance of an ideal array in this 
energy range. 

In case of the diffuse source sensitivity the number 
of telescopes is the found to be the most important fac¬ 
tor. Again the addition of MSTs of either type (SC or 
DC) would result in a considerable improvement com¬ 
pared to M25DC (similar to Array I) in the whole en¬ 
ergy range. However the improvement is slightly less 
significant when compared to the relative improvement 
in the point-source sensitivity. 


5. Conclusions 

This paper describes a new simulation and analy¬ 
sis chain that is used to study and compare array and 
telescope design concepts for CTA. We specifically fo¬ 
cus on the role of MST arrays which are optimized for 
performance in the core energy range of CTA between 
100 GeV and 1 TeV. The simplified detector model used 


for this study allows for investigation of a wide range 
of telescope parameters: effective light collection area, 
optical PSF, camera pixel size, effective camera trig¬ 
ger threshold, and effective integration window in time. 
The simplified telescope description allows us to isolate 
the most important telescope design characteristics and 
fully explore their influence on the performance of the 
full array. Realistic telescope designs can be mapped to 
our simplified detector model by choosing telescope pa¬ 
rameters that are matched to the physical characteristics 
of each design (mirror area, focal length, photosensor 
efficiency, etc.). This paper also investigates several as¬ 
pects of the array geometry optimization including the 
impact of the number of telescopes and their separation 
on array performance. 

A benchmark telescope array was used to assess the 
influence of each of the telescope and array parame¬ 
ters. Performance is evaluated for nominal observing 
conditions corresponding to a zenith angle of 20° and 
an NSB rate computed for a dark extragalatic field. We 
also examined the influence of the GL and higher NSB 
rates. Under all conditions, an optimized analysis is per¬ 
formed using a likelihood reconstruction based on sim¬ 
ulated image templates and BDTs for signal extraction. 

The likelihood reconstruction based on simulated 
templates offers a factor of two improvement in point 
source sensitivity (30-40% improvement in gamma-ray 
PSL), as well as a reduced energy threshold relative to 
image moment-based analysis. The likelihood recon¬ 
struction takes advantage of the possibility of fully re¬ 
solving showers with a finely pixelated camera. This 
technique, coupled with BDTs for event selection, al¬ 
lowed us to compare arrays very close to their maximum 
achievable sensitivity. 

We find that the substantial improvements in both the 
gamma-ray point-source sensitivity and angular reso¬ 
lution of an IACT array can be realized by telescopes 
with imaging resolution better than current-generation 
IACT designs. We find a 3CM-0% improvement in the 
gamma-ray point-source sensitivity between 100 GeV 
and 3 TeV when the telescope pixel size is reduced 
from 0.16° to 0.06°. The gain in point-source sensitivity 
comes primarily from the improvement in the gamma- 
ray angular reconstruction enabled by the higher reso¬ 
lution imaging of the shower axis. Over the same en¬ 
ergy range, the performance of an MST array is much 
less sensitive to the telescope light collection area and 
trigger threshold. We find that these parameters are im¬ 
portant in determining the array energy threshold but 
have little influence on the array performance above the 
threshold energy. 

With higher resolution shower images, the GL be- 
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comes more relevant than ever for the sensitivity of an 
IACT array. To determine the impact of the GF, we 
compared the same array simulated with values of B ± 
between 0 pT and 20.7 pT. For an MST array, the im¬ 
pact of the GF is largest around 100 GeV where the 
point-source sensitivity is reduced by 50%. The GF 
should be an important factor in selecting a site for 
future arrays and possibly for designing an observing 
strategy. 

Increasing the number of telescopes in the array ex¬ 
pands the effective area, improves reconstruction, and 
increases background rejection capabilities. The sensi¬ 
tivity can be improved faster at very low and very high 
energies by adding LSTs and SSTs. However, in the en¬ 
ergy range between a few hundred GeV and tens of TeV, 
expanding the MST array efficiently improves the sensi¬ 
tivity, regardless of the telescope design. In the limit of 
a finite array for which uncontained showers constitute a 
significant fraction of the total reconstructed event sam¬ 
ple, the improvement in point-source sensitivity scales 
faster than the square root of the number of telescopes 
between 300 GeV and 3 TeV. If the baseline CTA design 
is expanded to include 36 more MSTs, the point-source 
sensitivity in the core energy is improved by a factor of 
two. 

When considering arrays with the same number of 
telescopes, we find that the SC telescope design yields 
a 30-40% improvement in point-source sensitivity over 
the DC telescope design because of its superior imag¬ 
ing resolution. The DC telescope on the other hand 
has a slightly lower energy threshold resulting in better 
point source sensitivity below 75 GeV. The improved 
performance in a wide energy range from the SC design 
warrants further investigation. The improved sensitiv¬ 
ity reduces the total exposure time required for every 
science topic, while the smaller gamma-ray PSF addi¬ 
tionally helps with source confusion and morphology 
studies. The higher resolution shower images of the SC- 
like telescopes are also much less affected by noise from 
NSB. This translates to a much lower energy thresh¬ 
old during brighter sky conditions, e.g. in the galactic 
plane. This may lead to a much higher effective duty 
cycle since observations can be continued into brighter 
moon phases without sacrificing the low energy regime. 

While the SC-like array is more sensitive in compari¬ 
son to the DC-like design, no SC MST telescope has yet 
been built. Construction of an SC prototype at the site of 
VERITAS is under way. This prototype offers a chance 
to study the performance of the SC optics in realistic 
circumstances. This experience should also provide a 
more realistic cost model for the two-mirror systems. 

At this point in the design of CTA, it is unlikely that 


all MST telescopes would be of the SC design. If the SC 
prototype can be built successfully and cost-efficiently, 
the baseline CTA array could be expanded to include an 
additional number of SC MSTs. The study of mixed ar¬ 
rays is ongoing. No matter which optical design is cho¬ 
sen, expanding the MST arrays offers significant ben¬ 
efits for the performance of CTA in the central energy 
range between 100 GeV and 1 TeV. 
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